arXiv:1504.03288vl [hep-th] 13 Apr 2015 


Perturbatively Charged Holographic Disorder 


Daniel K. O’Keeffe0 and Amanda W. Peet@ 

Department of Physics, 

University of Toronto, 

Toronto, Ontario, 

Canada MSS 1A7. 


Abstract 

Within the framework of holography applied to condensed matter physics, we study 
a model of perturbatively charged disorder in D = 4 dimensions. Starting from initially 
uncharged AdS^, a randomly fluctuating boundary chemical potential is introduced by 
turning on a bulk gauge field parameterized by a disorder strength V and a characteris¬ 
tic scale ko- Accounting for gravitational backreaction, we construct an asymptotically 
AdS solution perturbatively in the disorder strength. The disorder averaged geome¬ 
try displays unphysical divergences in the deep interior. We explain how to remove 
these divergences and arrive at a well behaved solution. The disorder averaged DC 
conductivity is calculated and is found to contain a correction to the AdS result. The 
correction appears at second order in the disorder strength and scales inversely with ko. 
We discuss the extension to a system with a finite initial charge density. The disorder 
averaged DC conductivity may be calculated by adopting a technique developed for 
holographic lattices. 
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1 Introduction 

The AdS /CFT correspondence has proven to be a remarkably powerful tool for probing the 
detailed structure of strongly coupled quantum held theories. Its broad list of successes 
includes applications to modelling the quark-gluon plasma, condensed matter phenomena 
such as superconductivity, and even huid dynamics. Comprehensive reviews of these subjects 
include, PO. 0 and |3], respectively. 

An underlying theme to this progress is the reduction of symmetry in holographic mod¬ 
els. Systems which possess too many symmetries display behaviours that are not desirable 
in condensed matter models, a prime example being the inhnite DC conductivity dual to 
the Reissner-Nordstrom-AdS (RN-AdS) geometry. This is not a surprising feature. The un¬ 
derlying RN-AdS geometry possesses translational invariance along the boundary directions. 
This means that the charge carriers in the dual theory have no means by which to dissipate 
momentum, resulting in an inhnite DC conductivity. Realistic condensed matter systems do 
not display this behaviour, so if a holographic model is going to be useful for studying these 
kinds of problems, we need a way to break translational invariance. 

Several avenues of investigation have been carried out, including explicit holographic lat¬ 
tices: 111 , 0 , 0 , 0 , 0 , 0 , 1101 . 1111 , Q-lattices: [12] , [13] , [HI and breaking diffeomorphism 
invariance: |15j . ra. im. [T8] . [T9] . [20] . [2T] . [22] . Other examples of holographic symmetry 
breaking include the breaking of rotational symmetry m, m, E5], m, m, relativistic 
symmetry [2B], [29], [50] and hyperscaling violation [5T], [32], [35], [51]. Our focus will be on 
disorder. Within the context of bottom-up modelling, we will report here on a perturbative 
construction of holographic disorder. That is, we will seek a particular spacetime solution 
that is sourced by a disordered, randomly huctuating held in the dual theory and use it to 
study the dual DC conductivity. 
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Disorder is a common feature in real world condensed matter systems, but it is difficult 
to model using traditional field theory approaches: little is known at strong coupling. This 
is especially true in the context of localization. It is well known that in a non-interacting 
system, the addition of disorder can completely suppress conductivity [35]; for a compre¬ 
hensive review, see [36]. Turning on interactions complicates the situation and many-body 
localization may occur [57]; a review may be found in [55] • Experimental studies of many- 
body localization are also in their infancy. Since holography is a strong-weak equivalence, 
it provides a potential way forward for studying disorder. Understanding what the gravity 
dual tells us about localization might even provide a definition of many body localization. 
Exactly what this would look like in a gravitational theory and how the bulk fields would 
behave is an open and interesting question. 

Several approaches to disorder within holography have been proposed. See [55], [ID] and 
[TT] for early studies. The replica trick for disordered CFTs was extended to holography in 
[52] and later in [55] . 

To fully see the effect of quenched random disorder, [55] proposed that disorder should 
be modelled by applying random boundary conditions to a bulk held and allowing the held 
to backreact on the geometry. The source of disorder is characterized by a distribution 
Pi/[hE(a;)], over random functions W{x). The subscript V is in reference to another function 
V{x), which is taken to be the boundary value of a gauge held in the bulk. The system 
is assumed to be self-averaging, meaning that over length scales much larger than than a 
typical disorder scale, homogeneity is restored. An interesting feature of this construction 
is that the entire functional Py[lU(a;)] runs as the energy scale changes. To this end, ]55] 
and later ]55], worked to construct a holographic functional renormalization scheme so that 
disorder averaged thermodynamic quantities could be computed. 

Another approach to disorder in gauge/gravity duality was proposed in ]56]. A particular 
background which is already deformed away from AdS, but still satishes the null energy 
condition, is considered. The degree to which the geometry is deformed is controlled by a 
parameter which is meant to represent the amount of disorder in the system. The wave 
function of a probe scalar in this background is studied by looking at the nearest neighbour 
spacing distribution of the pole spectra of the two-point correlator. By increasing the amount 
of disorder, a transition is observed between an initial distribution and a Poisson distribution, 
which is analogous to what happens in a disorder driven metal-insulator transition. 

More recently, mi applied a spectral approach to modelling disorder in a holographic 
superconductor. The basic idea is similar to [55] where the disorder is sourced by a random 
space-dependent chemical potential by setting the boundary condition on a bulk f/(l) gauge 
held. In this proposal, the chemical potential takes the form of a sum over a spectrum which 
depends on one of the boundary directions, x 



( 1 . 1 ) 


k=ko 


where Sk is a function of the momenta k which controls the correlation function for /i(x); 
different choices of Sk lead to a different values of the disorder distribution average. U is a 
parameter which sets the strength of the disorder, ko and fc* dehne an IR and UV length 
scale cutoff for the disorder, respectively. 6k is a random phase for each value of k. A spectral 
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representation like this is known to simulate a stochastic process when fc* is large [IH]. 

The dirty chemical potential (11.11) is incorporated into an Einstein-Maxwell model. The 
Maxwell equations are solved numerically in an electric ansatz with the boundary condi¬ 
tion that At approaches fll.ip near the boundary. The analysis is done in a hxed AdS- 
Schwarzschild background. Evidence is found for an enhancement of the critical temperature 
of the superconductor for increasing disorder strength. Numerical evidence that localization 
occurs in this model was reported in [1^. This model has also been extended to holographic 
p-wave superconductors in [50] where the same behaviour of the critical temperature with 
disorder strength was observed. 

There is a resemblance between this spectral approach to disorder and the construction of 
a holographic ionic lattice in [B] and [BT] where the lattice is set up by a periodic boundary 
chemical potential. The difference is that for a holographic lattice, there is only a single 
periodic source of a hxed wavelength. In the spectral approach to disorder, there is a sum 
over periodic sources of arbitrary wavelength, so the effect of disorder may resonate more 
strongly throughout the entire bulk geometry, having non-trivial effects deep in the interior. 

Recently, [52] applied a spectral approach to modelling disorder sourced by a scalar held 
in 2-1-1 bulk dimensions. The initial clean geometry in this case is AdS^ and the scalar held is 
allowed to backreact on the geometry in the spirit of [H]. By treating the disorder strength 
as a perturbative handle, a second order analytic solution for the backreacted geometry is 
obtained. It is observed that the disorder averaged geometry, in the deep interior, takes the 
form of a Lifshitz metric with a dynamical critical exponent z set by the disorder strength. 
Numerical solutions are constructed for strong disorder, also display this behaviour. A 
similar implementation of scalar disorder was used in [53] and [M] to study the conductivity 
of holographic strange metals with weak quenched disorder. 

In [To], it was shown that a wide variety of inhomogeneous IR geometries arise as solutions 
to Einstein-Maxwell theory with a single periodic source at the boundary. These observations 
suggest that studying disordered boundary potentials in Einstein-Maxwell theory at both 
zero and hnite net charge density may lead to novel gravitational solutions and insight into 
disordered condensed matter systems. 

Our goal here is to study charged disorder with backreaction. As we will see, the general 
problem for a system with baseline charge density is very complicated. Accordingly, we 
study the story perturbatively. We build in the disorder the same way as m, ISO], with 
a randomly varying boundary chemical potential, and we also include backreaction. In the 
bulk, this corresponds to having a fluctuating gauge held and letting it backreact on the 
initially clean AdS geometry. We then solve the Einstein equations perturbatively in the 
strength of the disorder and construct an analytic, asymptotically AdS, solution. Similar to 
the case of scalar disorder studied in [52], the disorder averaged geometry contains unphysical 
secular terms which diverge in the deep interior. We explain how these divergences may be 
tamed and ultimately hnd a well behaved averaged geometry. 

Our primary interest is in investigating the transport properties dual to the disordered 
geometry. Adapting a technique hrst developed in [13], we directly calculate the disorder 
averaged DC conductivity. We hnd a correction to the conductivity dual to pure AdS starting 
at second order in the disorder strength which scales inversely with the smallest wavenumber 
ko in the disordered chemical potential. We will also discuss extensions to systems with hnite 
charge density. 
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Holographic momentum dissipation and disorder has also been approached with the con¬ 
text of massive gravity. See [55] , [56] , isa, [58] , [59] , [60] , and [6l] . While this approach will 
not be our primary focus here, it will nevertheless be interesting to understand how results 
within massive gravity mesh with explicit implementations of holographic disorder. 

This paper is organized as follows. In section |2] we construct a bottom-up, disordered, 
holographic spacetime analytically. Using the strength of the disorder as a perturbative 
handle, we solve the bulk equations of motion up to second order. Using a Gaussian distri¬ 
bution for the disorder source, we explicitly evaluate the disorder average of the backreacted 
solution. In section [3l we implement a resummation procedure for removing spurious secu¬ 
lar terms which develop in the disorder averaged geometry as a result of our perturbative 
expansion. With the regulated solution in hand, we study the resultant DC conductivity in 
section m Section [5] points out some features of incorporating disorder in a geometry with 
an initial charge density. We point out how techniques used to study transport properties 
of holographic lattices may be used to access the DC conductivity in the hnite charge case. 
Finally, in section [6] we summarize our hndings and comment on possible directions for future 
work. 


2 Perturbatively charged disorder 


We are interested in studying the effect of a disordered holographic lattice, with the ultimate 
goal of understanding transport properties like conductivity in the dual model. To this end, 
we start with a gauge held the D = d + 1 dimensional bulk and then introduce a random 
perturbation around the initially clean baseline solution. The model is then Einstein-Maxwell 
gravity 


5 = 


J + 


d{d — 1) 

L2 




( 2 , 1 ) 


IbvrGjv 

Note that since we want an asymptotically AdS solution, we have set the cosmological 
constant to A = —d{d — 1)/2L^, where L is the AdS radius. The equations of motion we 
need to solve are then the Einstein and Maxwell equations. In what follows, it will be 
convenient to work with the traced Einstein equations, that is we need to solve 



1j 



1 

2(d- 1) 



( 2 . 2 ) 


= 0. (2.3) 

We consider a system where the initial charge density is zero. In this case, the baseline 
solution is just AdSd+i with At = 0. That is, the initially clean (no disorder) background is 

ds^ = ^ (—-I- dr^ -I- dxf) , (2.4) 

where the boundary is at r = 0 in these coordinates. To introduce disorder, we will perturb 
about this background by turning on the gauge held. To get a tractable system, we will 
introduce disorder along only one of the boundary direction x. According to the holographic 
dictionary, near the boundary, the time component of the gauge held behaves like 


At{x,r) = iJ,{x) + p(x)r‘^ ^ H-, 


(2.5) 
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where /r(a;) is the chemical potential and p{x) is related to the charge density of the dnal 
theory. We sonrce the disorder by taking the chemical potential to be a snm of periodic 
fnnctions with random phases nsing a spectral representation 

N-l 

p{x) = V'^An COs{knX + On) ■ ( 2 . 6 ) 

n=l 


Here, H is a constant which controls the strength of the disorder which we will nse as a 
pertnrbative handle. Similar profiles for the gange field have been stndied in the context of 
disordered holographic supercondnctors in [17] , |19], [50]. A spectral representation for scalar 
disorder was used in [52]. The wavenumbers kn = nA/c are evenly spaced with Ak = k^/N. 
1/fco is thought of as a short distance cutoff for the disorder and is held fixed so that in the 
limit that A —)■ oo, Afc —)■ 0. The On are random angles that are assumed to be uniformly 
distributed in [0,27r]. The amplitudes are 

An = 2^S{kn)Ak, (2.7) 


where S{kn) controls the correlation functions. Since the On are uniformly distributed, aver¬ 
aging over them (the disorder average) is done simply by taking 


{f)D - 


l.m /n-/ 

N^+oo I 27r 
'' i=i 


( 2 . 8 ) 


If the function S{kn) is taken to be 1, this leads to p{x) describing Gaussian noise. In other 
words, by explicit calculation, we have that 

_2 

{p{x))d = 0, {p{xi)p{x2))d = V 6{xi-X2). (2.9) 


The gauge field At sources the disorder in the system, which means that it contributes at 
order V. This in turn implies that the right hand side of the traced Einstein equations fl2.2p 
start contributing at second order in the disorder strength. Only Maxwell’s equations receive 
any corrections at first order in the disorder strength, so we can solve for At in the clean 
AdS background. The solution which obeys the boundary condition fl2.6p at r = 0 is 


A Ad-2)/2 

At(r, x) = 1 / ^ cos(knX + On) , (2.10) 


where K(^d- 2 )/ 2 {knr) is the modified Bessel function of the second kind. Notice, by turning 
on disorder, we went from a solution that had zero charge density to one that has a finite 
charge density. We can compute this in the dual theory: it is (</*), i.e. the expectation value 
of the current density. The result is 


(J‘) = 


{d - 2)VL’^-^ 


SttG 


N 


^Ankt^m 

n=l 


d)/2) 


2 '^-! r((d- 2 )/ 2 ) 


COs{knX + Or, 


( 2 , 11 ) 


Note, this formula is valid when (2 — d)/2 is not a negative integer (which will be the case in 
what follows). Now, we can plug the solution for At back into the traced Einstein equations 
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and work out the metric corrections at second order in V. This is generally difficult to 
do in arbitrary dimensions. We will focus on the case with D = 4, in this way the dual 
to the initially clean geometry is a 2 + 1 dimensional CFT. We need an ansatz for the 
backreacted metric. Since we are adding a perturbation in the x direction only, the metric 
coefficients could generally depend on both the radial direction r and the boundary direction 
X. Furthermore, we insist that the geometry is asymptotically AdS and we will work within 
Fefferman-Graham gauge. At the end of the day, this means that the general form of the 
backreacted metric is 

ds'^ = — [—a{x,r)dt‘^ + dr'^ + r]{x,r)dx‘^ + S{x,r)dy'^) , (2.12) 

where the functions a{x,r), r]{x,r) and 6 {x,r) need to be solved for. We solve for them in 
a perturbative expansion 


_2 

a{x, r) = l-\-V 0 : 2 ( 3 ;, r) + • • • 

_2 

ri{x, r) = 1 + V 772 ( 3 ;, r) + •• • 

5(3;,r) = 1 + v‘^62{x,r) H-, 

so to second order we solve for 02(3;, r), 772(3;, r) and 62 {x,r). The traced Einstein equations 
(I2.2j) expanded to second order give the following set of equations for the metric coefficient 

- 2rdrS2 + 2r^dla2 - 2rdrTi2 + 2r^d'^a2 - Qrdra 2 = [{drHY + , (2.13) 

4 

-2r‘^dl52 + 2rdrr]2 + 2rdra2-2r‘^dla2-2r‘^dlr]2 + 2rdr52 = [{drllY - [d^HY] , (2.14) 

2 

5r4(«2 + ^ 2 ) = ^{drHYd^H ), (2.15) 

4 

Qrd 2 r ]2 + 2 rdra 2 + 2 rdrS 2 - 2r^cl^02 - 2 r‘^dlS 2 - 2 r^d‘^r ]2 = ^ [{drllY - [O^HY] , (2.16) 

4 

2rdrri2 + 2rdra2 + Qrdr 62 - 2r‘^d^62 - 2r‘^dl52 = ^ [{drliY + [O^HY] , (2-17) 

where H = H{x,r) such that = VH{x,r) is the solution to the Maxwell equations at 
hrst order fl2.10p . 
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2.1 Second order solution 


By using combinations of equations fl2.13p - fl2.17p . we can solve for a 2 {x,r), ri 2 {x,r) and 
62 {x,r). The solutions are 


«2(x,r) = -H^{x) + -Cl + -Car - + Wn + 2fcnr + 1)< 


, ‘ZkfiV 


2 
N-l 


2 2 


16L2 ^ kl 32L2 ^ kl 

n=l n=l "■ 


1 

+ 22 2knr 1) exp(-2fc„r) cos^{knX + 9n) 

n=l ^ 

1 4 4 

+ qT 2 XI X + 2 knkmr‘^] COs{{kn - km)x + 6 n- 6 m) 

SL k^km 

n=l m=l 


N-l n-l 

1 

Q r2 / j / j u u 
n=l m=l ^ ^ 


1 A ^—1 n—l A A 1 1 

1 X ^ X ^ -^n^mknkm 


EE 


n=l m=l 


[1 -f {kn - km)r] e cos((/c„ - /Cm)a; -h 6^^ - 6^^) 

[(^n + km)‘^T'^ + 2[kn + km)^ -|- 2] C 


X COs((/u 7 ^ ~\~ k^n^X X Ofi ^m) ; 


I 2^71, T* 


1 11 1 /12 S ^ ^ 42 

= -H,i.) + -c. + + _ 5 ; ^ _ _ 5 ; + 2 ,„. + 1 ). 

rt=l “ n=l 

1 ^42 

+ X exp(-2/c„r) cos^(/c„a; + 6'^) 

n=l ^ 

1 n—l . . 

^EE^ — [l -F (fcn -F km)r -F 2knkmr‘^] e COs{{kn - km)x + 9n- 9m) 


8L2 

n=i m=i 
AT-l n-l . . 

^ X X XX (,-(kn-km)r cos{{kn - km)x + 9n - 9m) 

oL/ 

n=l m=l 

1 A^~l n—l A A 1 1 

4ee#^ ^ [(^n + -|- 2(/c„ -|- km)T -|- 2] C 


2 ^^ n=l m=i 


X COs((/c^ “ 1 “ km)x “t“ 9^1 9 m) 


1 1 \ 

772 ( 0 :, r) = -r^a^i/dx) - Car^ + H^^x) - ^ + 1) exp(-2A;„r) 


4 

7V-1 


32L2 ^ A:2 

n=l 


1 "X-X 

X X(2^nr -h 1) exp(-2/c„r) cos^{knX + 9^) 


16L2 ^ kl 

n=l "- 

.. Af—1 n—l A A U U 
i X ^ X ^ ^n^m^n^m 




[(^n + km)r -F 1] exp(-(/c„ -F km)r) COs((/Cn + km)x + 9n + 9m) , 
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where ci and C 2 are constants and Hi{x) and H^lx) are arbitrary fnnctions of x. The solntion 
for the metric coefficient in the x direction, ri 2 {x, r), is the simplest, as the eqnations of motion 
fl2.13p - fl2.17p can be arranged to give a hrst order differential eqnation for ri 2 {x,r). 

We want a solntion that is asymptotically AdS as r —)■ 0 and is regnlar in the interior. 
The former condition is tantamonnt to insisting that a 2 {x, 0) = ri 2 {x, 0) = 62 {x, 0) so that the 
metric takes on the correct asymptotically AdS form in Fefferman-Graham gange. Imposing 
these conditions hxes C 2 = 0 and Hi{x) = const and H^lx) = const. The solntions for 
a 2 {x,r), r] 2 {x,r) and 62 {x,r) which are asymptotically AdS are regnlar in the interior are 


, ,1 I ^Al I ^Al 
2 ^ 16L2 32L2 ^ 

n=l “ n=l 


E T? + WI E 7?+ 2'=”’' + 1) <=>‘P(-2M.O 


1 A^ 

^ -^{2klr^ 2knr 1) exp(-2A;„r) cos^{knX + 9n) 


16L2 ^ kl 

71=1 

N I n-1 . . 


+ ^ V] V] -TT^ [l + {kn + km)r + 2knkmr‘^] e COs((fc„ - km)x + 0^ - 9m) 

' kjlkyyi 
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n=l m=l 
N—1 n—1 


~ ^ ^ COs((fc„ - km)x + 9n- 9m) 


71=1 771=1 

N—1 n—1 


h h 

rvn.rvr) 


^ E E [(*=« + + 2 (*'« + Ur + 2 ] 


2 T 2 (^kn + km)'^ 

71=1 m=l ^ ' 


X COs(^(^kfi “1“ km)x “t“ 9n “1“ 9m) 7 


3 ^'^2 


1 1 A'^ 

62 {x, r) = -c+ ^ TT(2fcy + 2knr + 1) exp(-2fc„r) 


2 16L2 ^ kl 32L2 ^ kl 

n=l n=l 


1 A'^ 

^ -j{2klr‘^ + 2knr 1) exp(-2/c„r) cos^{knX + 9n) 


16L2 ^ kl 

71=1 

N I n-1 . . 

X ^ X ^ ^71^771 


s^EE 


n=l 771=1 

N—1 n—1 


h h 

^71^771 


[l d- {kn + km)r + 2knkmr‘^] e COs((fc„ - km)x + 9n- 9m) 


^ ~ ^ COs((fc„ - km)x + 9n- 9m) 


8L2 ^ ^ k„,k 

71=1 m=l 


A A h Jc 

■^ 71 -^ 771 ^ 71^771 


^ A ^—1 n —1 

2T2 ^ ^ {kn + kmY 

71=1 771=1 ^ ' 


[(^n + + ‘^{kn + + 2 ] e 

X COs(^(^kfi km}x 0>fi ^rn) 7 




















1 

V 2 {x,r) = + 1) exp{-2knr) 

n=l 

1 ^'^2 

+ Z-/ exp(-2/c„r) cos^(/c„a; + 6 ^) 

n=l ^ 

1 n—1 A A 1 1 

+ 72 exp{-{kn + km)r) COs((/c„ + km)x + 6 n + Om) , 

where c is a constant. The fact an asymptotically AdS solution exists was to be expected 
as in d + 1 = 4 dimensions, the gauge held which encodes the disorder sources a relevant 
operator in the dual theory. Geometrically, the gauge held falls oh sufficiently fast as r 0 
so that it does not disrupt the asymptotic form of the metric. 


2.2 Disorder average 

The solutions we have found for a 2 {x,r), r] 2 {x,r) and 62 {x,r) simplify considerably under 
disorder averaging. Using fl2.8p . we average over the random angles 6 n- The resulting sums 
may be expressed in terms of special functions 


(q! 2 )d = lim 

A^^+CO 


2^ 4fcoL2 


+ 2fcorln(l - exp(-2fcor/iV)) 


(2.18) 


+2/co?"exp(—2fco^")d>(exp(—2^0?"/-^), 1, N) — NLi2 [exp{—2kQr/N)] 
2ky exp{-2kor{N - 1)/N) - 1 


N 


exp{2kor/N) — 1 


+ N exp{-2kor)^{exp{-2kor/N), 2, N) 


{ 62 ) 0 = lim <-c + 


:NH, 


( 2 ) 


jvGToo ^2' ' 4A;oL2^’^ Ak^L'^ L 

—2A;o?"exp(—2A;o^")‘h(exp(— 2 ^ 0 ?"/-^), 1,-^) + NLi 2 [exp(—2fco?"/^)] 
2klr‘^ exp{-2kor{N - 1)/N) - 1 


— 2/co?"ln(l — exp(—2A:o?"/-^)) (2.19) 


N 


exp(2fcor/iV) — 1 


— Nexp{—2kor)^{exp{—2kor/N), 2, N) 

{V 2 )d = c (2.20) 

where — is the generalized harmonic number, Li 2 {z) = — is the poly- 


n=l 


71=1 




logarithm of index 2 and <h(z, a, b) = ^ is the Lerch Phi function. By taking the 


^ (b + nY 

n=0 ^ ' 

large N limit, the above expressions are found to be hnite and give 

1 




(1 - V)exp(-2/cor), 


4A.„n 


(1 - V) + 77^(1 - V) exp(-2V), 


4fcoL2 


( 2 . 21 ) 

( 2 . 22 ) 
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{V2)d = c, (2.23) 

The disorder averaged solutions ( 02 )d and ( 52 )d are divergent in the interior of the geometry 
as r —>■ 00 . Since the gauge field which sources the disorder is a relevant perturbation, this 
is not surprising. In order to have a hnite solution in the interior, we should resum the 
perturbative solution in the spirit of the Poincare-Lindstedt method for removing secular 
terms (terms that grow without bound) in perturbative solutions to differential equations. 
This kind of situation was also encountered in the case of scalar disorder studied in [52] . 


3 Resummation of disordered solution 


We have seen that the disorder averaged metric coefficients ( 02 ) d 112.2ip and (^ 2 ) d 112.22p 
diverge in the deep interior as r —>■ cxo due to the presence of secular terms. This indicates 
a break down in the perturbative scheme and needs to be corrected. This outcome was not 
totally unexpected, as a similar sort of divergence in one of the disorder averaged metric 
components was found |52]. This divergence was corrected by resumming the disordered 
solution in the spirit of the Poincare-Lindstedt method. We will adapt this technique to our 
problem and see that a resummed solution, which has (q! 2 )d and { 62)0 finite in the deep 
interior, is available. 

The procedure is as follows. We look for additional terms that could contribute at the 

— 2 

correct order in the disorder strength, V , as they should only correct the unphysical secular 
terms in the disorder averaged metric. Furthermore, these new terms must not violate the 
asymptotically AdS condition imposed on the uncorrected disordered geometry. The simplest 
ansatz then for corrected backreacted metric in Fefferman-Graham coordinates is 






a{x, r) 


dt^ + dr"^ + ri{x, r)dx^ -f 


d{x,r) 

l32{rY(y) y 



(3.1) 


The new functions /9i(r) and /92(’") may be chosen to remove the secular terms in ( 02)0 and 
{^ 2 ) 0 - Hence, they should only contribute starting at second order. Their exponents, then, 
Wiy) and P{y) should be expanded in powers of the disorder strength as 


W{V) = W 2 V^ + + • • • , (3.2) 

p(V) = P 2 F' + py^ + • ■ ■ (3.3) 


where, as in the usual PoincarGLindstedt procedure, the Wi and Pi are constant coefficients 
that we will adjust in order to remove the secular terms in ( 02 ) 1 ) (I2.2ip and { 62)0 (I2.22p . In 
order to ensure that the spacetime geometry remains asymptotically AdS, we must enforce 
the boundary condition f3i{r) —)■ 1 and ^ir) —)■ 1 as r —)■ 0. Moreover, in order to ensure 
that the functions f3i{r) and f32{r) have a chance at removing the problematic terms in the 
disorder averaged metric components, we should have that /5i(r) and yy diverge as r — )■ cxd . 
This is the only way that a new term can compete with (and ultimately regulate) the already 
divergent secular terms. 

The placement of the new functions /3i{r) and ^ir) in the metric fl3.ip is not arbitrary. 
Expanded to second order in the disorder strength, the right hand side of the traced Einstein 
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equations fl2.2p is left unchanged. This is sensible since the source of the disorder, namely 
the gauge field (I2.inp is not being modified. Notice also that (3i{r) and ( 32 {r) are only 
functions of the radial coordinate r as their sole purpose is to correct the secular terms in 
the averaged metric components (q! 2 )d (I2.2ip and { 82)0 (I2.22p . There is no new function 
associated with the x direction as {ri 2 {x,r))D (I2.23P is already finite everywhere and does 
not require a correction. 

The off-diagonal [r, x] component of the traced Einstein equations fl2.15p is also un¬ 
changed. This is crucial, as this component of the traced Einstein equation allows us to 
solve for the combination a 2 {x,r) + 82 {x,r), and so the x dependence in this combination 
will also not be modified from the original solution. The only difference brought about by 
including [81 and (82 is the possibility of modifying the overall partial integration functions of 
the radial coordinate in 02 ( 2 :, r), 82 {x,r) and T] 2 {x,r). It is this modification that will allow 
us to remove the secular terms from the disorder averaged metric components. 

Finding the solution for the metric coefficients in fl3.ip to second order in the disorder 
strength goes through exactly the same way as for the original case. We look for an expansion 
of the form 


_ 2 

a{x, r) = 1 + V a 2 {x, r) + ■ ■ ■ 

_ 2 

ri{x, r) = 1 + V r] 2 {x^r) -\- ■ ■ ■ 

8 {x,r) = 1 -h v‘^ 82 {x,r) -I-. 

Plugging this expansion into the traced Einstein equations fl2.2p yields a set of partial dif¬ 
ferential equations 

6W 2 P 

- 2rdla2 + 28282 + 2drr]2 - 2rdla2 + 69 ^a !2 -- ^drf82 

Pi P2 

- [{drHy + {d^Hf] , ( 3 . 4 ) 


o 00 2 W 2 2 P 2 

2rd^r]2 - 25^02 - 2 dr 82 2rd^a2 + 2 rdl 82 - 2drr]2 5^/52 




PI 


PI 




P 2 


Pi 


Pi ' P2 

81(3^ = ^ [{ 8 rHf - { 8 ^Hf] , (3.5) 


8x8r{a2 + 82) 


L2 


{ 8 rH){ 8 ^H), 


(3.6) 


2Wo 


2 P 2 


28r82 - 2r8^82 - 2r8^r]2 28ra2 + 69^^72 - 2r8^a2 -- ^drP2 


Pi 


P2 


= 5 HdrHf - (d^Hf] , 


(3.7) 
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(3.8) 


6 P 21V 

2drT]2 - 2rdl62 + 29 ^ 0(2 - 2rd'^52 + Qdr62 - ^drP2 - ^drl^l 

^2 f3l 


+ ^-^81^2 - ^-^{dr^2f = ^ [{drHf + {d.Hf] . 


[32 


[31 


Notice, the [r, x\ component fl3.6p is unchanged compared to the original [r, x\ component 
fl2.15p . as promised. In fact, the new traced Einstein equations differ only from the original 
ones fl2.13p - fl2.17p by the addition of terms with derivatives of /9i and (^ 2 - In particular, 
this means that the x dependence of the metric coefficient 0 : 2 ( 2 :, r), ri 2 {x,r) and 62 {x,r) will 
be unchanged. Again, this had to be the case since the source for the disorder, namely the 
gauge field fl2.1Up has not be modified, nor have the right hand sides of the traced Einstein 
equations fl2.2p . 

The solutions for 02 ( 2 :, r), r] 2 {x,r) and S 2 {x,r) are found to satisfy 


1 A 2 

02 ( 2 :,r) + 52 ( 2 :,r) = Gi( 2 :) + G 2 (r) + — X] (l + 2 kr,r + 2 klr‘^) e 


‘2knT 


8L2 ^ kl 

n=l 


1 n—1 A A 1 1 

72 


\{kn + km)^'f'‘^ + 2r[kn + km) + 2 ] 


X COs((/Cn + km)x + 9^ + , 

where Gi{x) is an arbitrary function of x and 


1 ^Al 


G2ir) = G'2(r) - ^ ^ (l + 2A:„r + 2ky) 


—2fc„r 


16L2 ^ kl 

n=l "■ 


(3,9) 


(3.10) 


The condition (13.91) along with the partial integration function G 2 {t) (I3.10p are also true for 
the original disordered metric solutions up to a new, potentially different, function G 2 {r). It 
satishes 


rIT; 

Ti 

21^2 


rP. 

T 2 

rW2 


rdlG2{r) - 2drG2{r) = + ^9^/92 - ^5r/92 

P2 


rP2, 


^r[3l - ^idr[3i) - -^{dr[32) 




(3.11) 


This is the new piece that will allow us to tame the secular terms in the disorder averaged 
solutions. 

The equations of motion fl3.4l) - 03.81) may also be combined to get 


1 


02 ( 2 :, r) - 52 ( 2 :, r) = Giir) - ^ + 


1 ^'A2 


8 L2 ^ kl ' 8L2 (^ + + ‘^^y) 

n=l n=l 


, ‘2knT 


1 A^—1 n—1 A A 

1 AnAm 


4^2 ^ ^ 2knkmy cos((/i:„ - km)x + 9n- 9m)e 

n=l m=l ^ ^ 


iV-l n-1 . . 

X X ^ X ^ 


-4Ij7E 


71=1 m=l 


h h 

•^n.n/n 


[1 + {kn + km)r\ cos{{kn - km)x + 9^ - 9m)e , 


(3.12) 
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where the function G^lr) plays an analogous role to G 2 {r) in (I3.9p . it satishes 


rd'^G^^r) — 2drG4{r) = 

Finally, the solution for ri 2 {x,r) is 
1 


2hF2 
■ /Si 


2P, 


rPFo 


drf^l + -5-9r/S2- 


P 2 


/3f 


(3.13) 


1 

7]2{x, r) = -r^d^Giix) + G^ix) - Ge{r) - ^ (2fc„r + 1) 


, ‘ZkjiT 


32L2 ^ kl 

n=l 


1 n—1 . . 

EE -nV^r 


L2 


i {kn A km) 

1=1 m=l ' ' 


[{kn + km)r + 1] cos((fc„ + km)x + 9 n + dm)e 


1 ^'^2 


+ ^ -j^ {‘^knT + 1) cos^(fc„a: + 9n)e 

n=l ”■ 


—2k„r 


where 


Ge{r) = dr 


1 02 /=< , , ^P‘i(^ol2 ^P‘2f,2o rW2 2 o 

;r9.G2 + ^(a,ft) + 


(3.14) 


(3.15) 


Using these solutions, we will be able to write down a resummed backreacted metric, devoid 
of secular terms under disorder averaging. 


3.1 Resummed disorder average 

With the solutions fl3.9p . fl3.12p and fl3.14p at hand, we can impose asymptotically AdS 
boundary conditions and require that the secular terms in the disorder average vanish. To 
do this, we need a choice for (3i{r) and /32(^)- A natural choice is 


Mr) = ft ( r ) = exp , 


(3.16) 


with hF 2 = 1 = —P 2 - V has units Length^'^^, so to second order and are 

dimensionless. With this choice, the partial integration functions G 2 {r) (13.101) and G 4 (r) 
fl3.13p simplify to 


G2{r) 


di + -d2r^ , 


(3.17) 


and 


G4(r) 


1 -1 

d-s + 2*^4’" + 


1 

- V 

2L2 ’ 


(3.18) 


where di, d 2 , ds and d^ are constants. 

Since G 2 {r) appears in 0/2 + S 2 and G 4 (r) appears in 02 — ^ 2 , the two metric coefficients 
will be modified from their original values by 


-{G2 + G4) — -{di + ds) + -(d2 + dijr^ + 
2 2b 


(3.19) 
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in a 2 {x, r), while 


\(Gi - G 4 ) = \(d, - d,) + - dty - 


(3.20) 


appears in 62 {x,r). In 772 , Gq (I3.15p becomes simply 


Gfi{r) = dr 


-rdlG2{r) 




(3.21) 


Setting the constant ^2 = 0 recovers the original solution for ri 2 {x,r). As such, ( 772 )^ is ex¬ 
actly the same as in fl2.23p . This is precisely as expected, ( 772 )d was already finite everywhere 
and requires no correction. 

Finally, setting ds = ^4 = 0 the corrections to ( 012 ) 7 ? and { 62)0 are precisely those needed 
to remove the secular terms. We hnally have 

+ 47)75 - exp(-2i„r) , (3,22) 

= T ■ 47)75 + 47 )l 5 <^ ■ exp(- 27 „r). ( 3 . 23 ) 

and ( 772 )D is given in fl2.23p . These results constitute the corrected disorder averaged metric 
coefficients, all of which are now hnite everywhere in the bulk. 

Note also that the choice for IF 2 = —P 2 ensures that there are no curvature divergences 
anywhere in the bulk. Calculating K = through order V in the disorder 

averaged metric gives K = const. 

With the resummed solution, we can ask about transport properties of the disordered 
geometry. In particular, we are interested in understanding how the conductivity is modified 
from the initially clean AdS case. We tackle this question in the next section. 


4 Conductivity 

In order the calculate the conductivity, we need to turn on a perturbation to the gauge held. 
We will be interested in computing the conductivity along the disordered direction x. In pure 
AdS, the zero-momentum conductivity can be computed by turning on A^irR) = ax{r)e~^‘^^ 
and computing the retarded holographic Green’s function via the linearized bulk equations 
of motion. Taking the a; ^ 0 limit gives the DC conductivity which turns out to be a 
constant. This result persists even at hnite temperature and is due to bulk electric-magnetic 
duality [62]. The calculation is simplihed in the AdS case, since the gauge held perturbation 
does not contribute the to linearized energy-momentum tensor and so cannot source any 
new metric perturbations. 

The situation is not so simple when a non-zero background value for another component of 
the gauge held is turned on, for example in the RN-AdS geometry. This means that, generally 
speaking, the perturbation of the gauge held needed to measure the conductivity will couple 
to background gauge held component and contribute at linear order in the equations of 
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motion. As a consequence, new metric perturbations are sourced and can lead to complicated 
linearized equations. 

In the context of holographic lattices, such as those studied in 0.111 and [8] the problem 
is magnihed. Virtually everything that can be sourced is and the perturbation equations turn 
into a complicated set of partial differential equations which must be solved using numerical 
techniques. In [H], a holographic lattice, sourced by a periodic scalar held was studied 
analytically in a perturbative expansion about the lattice strength. Due to the perturbative 
nature of the lattice, the linearized equations of motion for the gauge held huctuation turn 
out to simplify considerably and only a single metric and scalar perturbation are sourced 
making the system more amenable to analysis. 

For our disordered case, the situation does not simplify quite so drastically, even though 
the disorder strength is assumed to be weak. The gauge held perturbation will mix with 
the background At, producing terms which are second order in perturbations (i.e. one power 
of the disorder strength V and Aj,). These will in turn source metric perturbations of the 
same order and in principle may source metric perturbations along every direction. We can, 
nevertheless, still access information about the DC conductivity of the system by employing 
a technique hrst proposed in [I3]. The idea is to turn on a perturbation that is linear in 
time. An analysis in linear response theory implies that the DC conductivity may be directly 
calculated by taking advantage of conserved quantities in the bulk. By using the equations of 
motion to relate the boundary current to the magnitude of the applied electric held, the DC 
conductivity follows from Ohm’s law. In [IT], this technique was adopted to study the DC 
conductivity of inhomogeneous holographic lattices at hnite temperature. In section [5], we 
will discuss how this technique may be applied to holographic disorder with non-vanishing 
initial charge density. 

We consider a perturbation to the gauge held A^ = a^ir) — Et, where E is the constant 
magnitude of the electric held in the x direction. At the level of the linearized equations of 
motion, this perturbation may further source metric perturbations {Sgu, Sgxx, ^9yyy ^9xt}- 
We have already elected to work in Feherman-Graham gauge when constructing the backre- 
acted spacetime in order to ensure an asymptotically AdS solution. As such, we will impose 
that the metric perturbations obey the usual Feherman-Graham expansion near the bound¬ 
ary. In other words, as r —)■ 0, 6gab = /'f'‘^)59ab{o){.x) + Sgab(i){,x) -|- C>(r^). Note that the 

metric perturbations will be at most functions of x and r and will be composed of combina¬ 
tions of periodic functions in x. This is because the metric perturbations are being sourced 
by terms made up of the gauge held perturbation and the original time component of the 
gauge held that sources the disorder, so there is no possibility to source any other kind of 
dependence. 

Using the bulk Maxwell equations V^E^'^ = 0, so in particular then the x and r compo¬ 
nents require that 

a. = 0 . (4.1) 

Hence, the current = ^J—gE^'^ is a conserved quantity in the bulk. As r —)■ 0, this dehnes 
the boundary current in the x direction. 

The conductivity can be read oh from Ohm’s law, namely = aE. Therefore, if we 
can hnd an expression for /E and take the disorder average of the result, we will have 
found the DG conductivity directly (up to an overall normalization constant set by the action 
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dJ). 

Our strategy will be as follows. First, we write down an expression for linearized about 

the gauge field and metric perturbations. Since this quantity is a constant, we can evaluate 

it anywhere in the bulk. A particularly convenient choice is at the boundary r = 0, where, 

in keeping with our choice of Fefferman-Graham gauge, we enforce asymptotic AdS falloffs 

for all of the perturbations. We then take advantage of our perturbative handle, the disorder 

strength V. The linearized equations of motion generally contain terms at many orders in V. 

We start with the leading order Maxwell equations which we solve for ax{r). This solntion 

is then fed into the next order in pertnrbations which conples ax{r) to order V terms. These 

source the metric perturbations, which may be solved for via the Einstein eqnations and 

expressed in terms of the magnitude of the applied electric field E. The results solve for a 

up to second order, after taking the disorder average. The linearized eqnations of motion 

_2 

contain higher order terms as well, snch as V Sgab- Since the metric flnctnations are already 
sonrced by terms of the form Vqx, these terms are already fourth order in pertnrbations and 
not relevant. 

The linearized expression for J®, expanded near the bonndary at r = 0 is 


r= [«(o)a,(i) + VH(^,)Sgtxio)] , (4.2) 

y a(o)h(o) 

where the gauge held perturbation has been expanded near the bonndary ax{r) = ax(o) + 
+ • • • • Also 


a{x,r) = 0 ( 0 )(a;) +ra{i){x) ^ - 

g{x, r) = r]{Q) {x) + rg^i) (x) H- 

5{x,r) = 5 ( 0 ) (x) +r5(i)(a;) H-, 

near r = 0. The fnnction H{x,r) is related to the disorder sonrce as At{x,r) = VH{x,r) 
and it has also been expanded near r = 0 as 

H{x, r) = H(^o){x) + rH(^i){x) H-. (4.3) 

The next step is to get an expression for the gauge perturbation ax{r). To do this, we nse 
the linearized Maxwell eqnations to leading order in pertnrbations. This is just the Maxwell 
equation in AdS, so the solution is easy to hnd 


ax^r) = hi + h2r , (4.4) 

Where hi and 62 are constants. This solution then conples to the disorder sonrce At and 

the original backgronnd metric coefficients to prodnce the metric hnctuations. In principle, 

the gange held hnctnation itself will then receive corrections from the metric perturbations, 

bnt this reqnires going to, at a minimum, third order in pertnrbation theory. The upshot is 

that we can use the baseline solntion for to extract the DC conductivity to second order 

in perturbations. The constant 62 can be hxed in terms of the applied electric held E. To 

_2 

see how, replace the time coordinate t with the ingoing coordinate v = t — r + 0{V ). Then 
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the full gauge field perturbation = ax — Et is a (finite) ingoing solution provided that 
62 = -E up to second order. This is exactly the condition required in pure AdS to recover 
the correct DC conductivity and it will ensure that our result contains the pure AdS case 
plus a possible correction, so that the hnal disorder averaged solution will take the form 

_2 

(o’)/) = CTAdS + (o'Disorder)n) • (4.5) 

That is, the original constant AdS DC conductivity plus a correction. 

With in fl4.2p we must hnd an expression for 6gtx{o) in order to get the conductivity. 
This can be accomplished via the linearized equations of motion. Despite being linearized, 
the equations of motion are still quite complex and finding a solution is a difficult task. This 
difficulty may be circumvented by making use of of perturbative expansion in the disorder 
strength. To get an expression for Sgtx{o) it suffices to solve the linearized Einstein and 
Maxwell equations up to second order in perturbations. That is, we keep terms of order Qx, 
V , Vqx and 6g. 

We have already solved the linearized Maxwell equations through second order in fl4.4p . 
No further corrections to the Maxwell equations are sourced until at least third order in per¬ 
turbations. The linearized Einstein equations work out similarly. There are no contributions 
to the linearized equations at hrst order in perturbations. This had to be the case as all of 
the metric coefficients and metric perturbations are at least second order perturbative terms. 
Every diagonal component of the linearized Einstein equations is just the original equation 
plus a correction due to a term proportional to VE which, in principle could source met¬ 
ric perturbations along these directions. Remarkably, none of the diagonal components of 
the linearized Einstein equations, nor the [x, r] component contain the metric perturbation 
6gtx that we need to compute the conductivity. In fact, the relevant metric perturbation 
decouples completely from the others and only shows up in [t, r] and [t, x\ components. The 
equations are 

rdxdrSgtx + 2rdr5gtx + Vr{dxH)E = 0 , (4.6) 

25gtx - r^dldgtx - 2rdr5gtx - Vr\drH)E = 0 , (4.7) 


with a solution 


VE s—y A VE 

6gtx{x,r) = — X] (2 + 2/c„r -F cos{knX + 9^) exp(-/c„r) -F 2—K{x) , (4.8) 

n=l 

where K{x) is a function of x that will be hxed momentarily. As promised, the metric 
fluctuation is sourced by a term that goes like V times a gauge held huctuation. With fl4.8p . 
we can pick off 5gtx(o) and get and expression for in terms of E. We now have all of 
the ingredients we need to get to the DC conductivity. Inserting the results back into fl4.2p . 
expanding to second order in the disorder strength 


^ = 1 + ^(«2(0) + 52(0) - h2{0)) - ^2 / . 

n=l ^ 


2V^ ^ 42 

-JY COS^ {knX + 9r, 


(4.9) 


2V^ ^ ! A A k 

COs{knX + 9n) COs{kmX -f 6m) + 2 — E[(i)K{x) , 


L2 


n,m=l 

m^m 


U3 
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where, as r —)■ 0 


a{x, r) = 1 + V‘^a 2 {x, r) -)■ 1 + ( 02 ( 0 ) + ^ 2 ( 1 ) H-) , 

r]{x, r) = 1 + V 772 ( 0 ;, r) ^1 + V ( 772 ( 0 ) + r 772 (i) H-) , 

6 {x, r) = 1 + V 52 {x,r) ^ 1 + V ( 52 (o) + ^"<^ 2 ( 1 ) + ■ ■ ■) • 

From the metric coefficient 02 , 772 and 62 we get 

0^2(0) + ^ 2 ( 0 ) - h 2 ( 0 ) = COS^iknX + 0^) - ^ 

n=l ^ n=l ^ 

1 n—1 A A U h 

+ E E (r;!”): <=“(('=»++»"+»-) ■ 

n=l m=l ^ ^ 

The last step needed to get the DC conductivity is to take the disorder average fl2.8p of 
fl4.9F To do this, we need to know the function K{x) appearing in fl4.8l) . This function may 
be fixed by requiring a finite large N limit for the conductivity as well as a finite value for 
^ 9 tx{o) and a vanishing disorder average for 6 gtx- The condition on the disorder average of 
5gtx actually follows from the requirement of a hnite large N limit for the conductivity. If 

{Sgtx)D was a constant as r —)■ 0, then the H(^i^5gtx{Q) term in fl4.2p would vanish under 

disorder averaging and would not be able to regulate the large N limit. The coefficient of 
K{x) can then be hxed by dimensional analysis. A good choice is 


K{x) 


2 ^-1 ^3 

^ -j COs{knX + On) , 

1 ^T) 

n=l 



(4.10) 


where ^( 77 .) is the Riemann zeta function. There is more than one choice for the coefficient of 
K{x), however every choice which is dimensionally consistent and satisfies our basic require¬ 
ments for 6gtx returns the same result for the disorder averaged conductivity. Applying fl2.8p 
to fl4.9p and restoring the units (which is due to the overall normalization we have chosen 
for the action fl2.ip □) we get 

1 

“ IGttGn ^ AnGNL^ko ' 

The disorder averaged conductivity takes the form (14.511 as claimed. The overall contribution 
due to disorder is a constant which comes in a second order in the disorder strength. This 
is sensible as the disorder induced corrections to the geometry (13.ip are at second order. 
In fact, our result echoes that of the single holographic lattice constructed in [51]. From 
the point of view of the conductivity, only the metric coefficients a{x,r), S(x,r) and ri{x,r) 
really matter. Under disorder averaging, these functions approach a constant towards the 
boundary r —)■ 0 which scales as U /(4/coT^), which is a natural dimensionless constant made 

^In pure AdS, the dual DC conductivity is anc = 1/e^, where the bulk gauge field is normalized by 
— l/4e^. In our conventions, e has been absorbed into the gauge field and we have pulled out an overall 
factor of 16ttGn. 
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up of the three scales in the system. It is therefore not surprising that the correction induced 
by the disorder source on the DC conductivity involves precisely this combination. Moreover, 
notice that as D —)■ 0, the average DC conductivity fid.lip reduces to the expected pure AdS 
result. 

Note that the result fld.lip is specihc to a disordered source with a Gaussian distribution 
fl2.9p . It is an interesting question as to whether or not changing the disorder distribution 
will affect the hnal result. In particular, the fact that the correction to the conductivity 
is positive and not negative as might be expected is a result of having chosen a Gaussian 
distribution. It is unclear that this will be the case if the disorder distribution is modified. 
Moreover, it is also unclear that a positive correction persists to higher orders in the disorder 
strength. It might be that further corrections come in with the opposite sign, but are only 
visible at larger disorder strength. 

It is tempting to try and address the question of using a different distribution starting from 
fl4.9p . However, to arrive at this formula, we have explicitly made use of the resummed metric 
fl3.ip which implicitly assumes a Gaussian distribution for the disorder source. In particular, 
the functions /9i(r) and /32(^) are specihc for this realization of disorder. Understanding how 
the form of the metric coefficients in the backreacted geometry is dependent on how disorder 
is implemented in the system and how this affects the dual transport properties are questions 
we leave for future work. 

5 Finite charge density 

The observations of the previous section all follow from a system with zero initial charge 
density. Turning on a non-zero background charge density can change the game drastically. 
The baseline, clean, geometry is now the Reissner-Nordstrom-AdS (RN-AdS) solution. In the 
RN-AdS background, the boundary directions preserve translational invariance, and so there 
is no way to dissipate momentum. Goupled with the fact that the solution no longer possess 
particle-hole symmetry like the pure AdS case, the net result is a divergent DG conductivity. 
From the perspective of our previous results in pure AdS, it is a natural question to ask how 
much we can say about adding disorder to this background. 

The action and equations of motion are the same as in fl2.ip . fl2.2p and fl2.3p . In this case, 
the initially clean geometry is Reissner-Nordstrom-AdS and the gauge held is not constant 
everywhere. The baseline solution is 



(5.1) 


There are d—1 boundary directions. The boundary is located at r = 0 in these coordinates, 
/(r) is given by 



{d - 2)/igrg 
2(d- 1)L2 



(d-2)^ 

2(d- 1)L2 



2 (d-l} 


(5.2) 



(5.3) 
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where tq is the location of the outer horizon and /io is a constant associated to the chemical 
potential of the dual theory [63]. The temperature is 


1 / (d-2)Vgrg \ 

dvrro V 2{d-l)L^ J ' 


(5.4) 


We introduce disorder the same way we did before by adding a perturbation to the gauge 
held and insisting on the boundary condition fl 2 . 6 p and regularity in the interior. As an 
example, if d + 1 = 5 and we take the limit where all the horizons coincide at r = rg (i.e. 
T = 0). The solution to the Maxwell equations at hrst order is 


At{x,r) = /io 1 - 


A 

+ 1 / 


r(i + %-)r(%-) 


n=l 


X 2F1 


r(A, 

A, 


1 + 2 - 


(rQ — r" 


2\ "I 


(5.5) 


A„ ^ L\n A 

— ,1 + —;A„; 
2 2 


3r2 


COs{knX + 6n) , 


/ _ 

where An = 1 + -w 1 H—Note that the solution (15.51) vanishes at the horizon r = ro as 
required for the gauge held to be well dehned. 

In d + 1 = 4, the solution is more complicated. 


At{x, r) = /io(l - )+^yi -^(^0 - 

roj ^ ^0 


iy/(2)(ro - r) - 2(ro + 2r) 


l-A„/2 


(5.6) 


X 



^n. 

1 . l + ^TT, 

2 i\/2{r—ro) 

2-^1 

2 ’ 

2 ’ 2 

’ i\/2{ro—r)—2{ro+2r) 

to 

A„ 1. 

2 ’ 2’ 

1-|-A„ . 2 iV 2 

2 ’ ^^/2-2 



/ _ 

where A„ = 1 + w 1 + . The solution fl5.6p vanishes at the horizon r = vq and obeys the 

correct boundary condition fl 2 . 6 p at r = 0 . 

The idea now is to plug this back into the equations of motion and work out the back 
reacted geometry to next order in V perturbations. This is more complicated than the 
initially uncharged case. Since the initially clean geometry requires a nontrivial gauge held 
from the get go, there is a baseline term in At{x,r) that survives on the right hand side of 
the equations of motion. In general this will mix perturbative and baseline terms, meaning 
that the equations of motion must be consistently solved at both hrst and second order in V. 
Also, the solution presented here is not in the Feherman-Graham coordinates of the previous 
section. This adds the complication that, in general, the backreacted geometry could have 
oh-diagonal Qra contributions. Since the perturbation is only along the x direction. At only 
mixes the r and x directions. Hence, the only oh-diagonal term that can be sourced is Qrx- 
This means that the back-reacted metric should take the form 


ds^ = 


dv^ 

a(x, r)dt^ + ^ + r]{x, r)dx^ + h(x, r)dy^ + 2 x(a;, r)dxdr 


V 


(5.7) 
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The perturbative expansion for the unknown functions this time is then 

a{x,r) = f{z) + Vai{x,r) + v‘^a 2 {x,r) H-, 

F(r) = /(r)+FFi(r)+FV 2 (r) + -- - , 

6{x, r) = 1 + V8i{x, r) + v‘^52{x, r) H-, 

r][x, r) = 1 + Vrii{x, r) + v‘^ri 2 {x, r) + • • • , 

X{x,r) = Vxi{x,r) + V^X 2 (x,r) H-, 

note that the expansion for x(x, r) starts at order V. This off-diagonal component is not 
present in the initial clean geometry, so in the limit that Id —)■ 0, this metric coefficient must 
vanish. Following the same logic as in section HI we should try and solve the traced Einstein 

equations up to second order in the disorder strength. At zeroth order, the traced Einstein 

equations are just those of the baseline RN-AdS solution, ensuring that the full solution will 
be expressed as the RN-AdS geometry plus corrections. This observation is what constrains 
the backreacted metric component grr to be only a function of r. The disorder is turned on 
along the boundary directions, so we do not expect it to modify the location of the horizon. 
At higher order, the traced Einstein equations are complicated and it is unclear if a compact 
analytic solution can be found. While we will not attempt to answer this question here, it 
would be curious to understand if the numerical techniques applied to holographic lattices 
and scalar disorder could be applied here. 

In lieu of a full analytic solution, we can nevertheless extract the form of the DC con¬ 
ductivity by applying the technique used for inhomogeneous holographic lattices in m- We 
will briefly review the salient features of this technique as we go along. We will make a 
modest assumption about the behaviour of the disordered solution near a horizon and see 
that, given the metric coefficients and gauge held solution, it is possible to extract the form 
of the disorder averaged DC conductivity directly in terms of horizon data. 

We will assume the form of the backreacted metric is (15.7p with the event horizon at 
r = ro and asymptotes to AdS^ near the boundary r = 0. As R —?• 0, the geometry must 
reduce to the RN-AdS solution. We will work at hnite temperature T and express the gauge 
held solution as 

At{x, r) = /io VH{x, r ), (5.8) 

where near the boundary, H{x,r) respects the disordered boundary condition (12.61) . In 
general, it will be a sum over periodic functions in x of arbitrary wavelength The same 
must be true of the metric coefficients in (15.Th . Note that, even though there are many 
periodicities in the disorder source, and hence the spacetime solution, there is a common 
periodicity 27rN/kQ. The metric coefficients and gauge held will also be functions of the 
random angles 9n in the disorder source (12.6|) . In order for the gauge held to be well dehned, 
it must vanish at the horizon, as in the RN-AdS case. Hence, the solution (15.Sh should vanish 
as r ^ ro. 

We will insist on an asymptotically AdS solution, so the metric coefficients in the back- 
reacted geometry (15.71) will need to fall oh appropriately near the boundary at r = 0. We 
will also need to insist on regularity of the solution at the event horizon, r = tq. The metric 
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coefficients in (15.71) shonld then be expanded near the horizon as 

a{x,r) =-47TT{r - To) [a(^o){x) + 0{{r - ro))-\ -] , (5.9) 

F(r) = -47rT(r - tq) [F(o) + 0{{r - ro)) H-] H-, (5.10) 

6{x,r) = 6^o)ix) + (r - ro)5(i){x) + C>((r - tq)^) + •• • , (5.11) 

r]{x,r) = V(o)ix) + {r - ro)v(i)ix) + 0{{r - rof) ^ -, (5.12) 

X(x,r) = X{o)(x) + {r - ro)x(i){x) + C>((r - ro)^) • • • , (5.13) 

H{x,r) = (r - ro)Hi{x) + 0{{r - tq)^) H-, (5-14) 

where each of the fnnctions in the near horizon expansion may themselves be expanded in the 


disorder strength V. For example a(o) = ao{o) + idai(o)(x) + v‘^a 2 (Q){x) + • • •, and similarly 
for the rest of (15.9p - (I5.14p . Note that in the limit that T —)• 0, there will be terms in a(x, r) 
and F{r) which are proportional to T~^ and stay non-zero. These terms originate from the 
backgronnd charge density in the initial RN-AdS geometry. Onr convention will be that 
the hrst snbscript indicates the power of V that mnltiplies the coefficients while the second 
snbscript in brackets indicates the order in the near horizon expansion. Note also that the 
fnnction H{x,r) in fl5.14p . which appears in the solntion to the gange held to hrst order in 
the disorder strength, only starts at O(r — ro). This is to ensnre that the gange held vanishes 
at the horizon, as reqnired by regularity. To keep the discussion as general as possible, as in 
m, we turn on every possible perturbation to both the gauge held and the metric, namely 
{at{x,r) ,ar{x,r), ax{x,r)} and {6gtt ,Sgrr ,Sgxx ,Sgyy ,6gj;r ,Sgtr ,Sgtx}- Note that we do not 
need to turn on a. y component to the gauge held as we will be interested in measuring the 
conductivity along the disordered x direction. To this end, we use a linear perturbation 

Ax{x,r) = ax{x,r) - Et, (5.15) 

where E is the constant magnitude of the applied electric held at the boundary. 

The perturbations will need to fall oh appropriately near the boundary so as not to 
destroy the AdS asymptotics. Moreover, we will insist on regularity at the horizon r = tq. 
This condition may be enforced on the perturbations by replacing the time coordinate with 
the ingoing Eddington-Finkelstein like coordinate 

V = t — (47rT)“^ ln(ro — r) • . (5.16) 

In light of the regularity condition, the constraints on the near horizon behaviour of the 
gauge held are 


at{x, r) = 0 ^( 0 ) (x) 0{r - Tq) ^ -, 

ar{x,r) = ^—{ar{g){x) + 0(r-ro)^ -) , 

fir) 

ax(x, r) = ln(ro - r) {ax(o) + 0{r - tq) ^ -) , 

along with the relations 


®r(0) 

®a:(0) 


E 

4 ^’ 


(5.17) 

(5.18) 

(5.19) 

(5.20) 

(5.21) 
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each of which follows as a consequence of regularity after switching to the ingoing coordinate 
fl5.16p . Whereas for the metric perturbations, regularity requires 


^9tt{x, r) = ^fir) {6gtt(o) (x) + • • •) , (5.22) 

^0 

Sgrr{x,r) = ^ {Sgrr{ 0 ){x) H-) , (5.23) 

^ 0 /V J 

5gxx{x,r) = Sgxx{o){x) + 0{r - ro) H-, (5.24) 

5gyy{x, r) = 55'ra{o) (^) + C’(r - ro) H-, (5.25) 

5gtr{x,r) = 5gtr{o){x) = 0{r - vq) ^ -, (5.26) 

5gxr{x, r) = (^fi'xr(o)(a:) + 0{r - ro) H-) , (5.27) 

5gtx{x, r) = 6gtx{o){x) + 0{r - ro) H-, (5.28) 

along with the conditions 

^9tx{o) ^9xr{o) 5 (5.29) 

^9tt{0) 4“ 9rr(0) ‘^9tr(0) 0 • (5.30) 


Notice that since the disorder source is a sum over periodic functions in x of arbitrary 
wavelength, the gauge held and metric perturbations will have to be as well. Again, they 
will all share the same common periodicity 27iN/kQ. Furthermore, the perturbations will all 
be functions of the random angles in the disorder source fl2.6l) . 

The next step, as in the case of a single holographic lattice, is to identify two useful con¬ 
served quantities in our background. The hrst is the boundary current, = yJ—gF^'' which, 
just as in section 0] is a bulk constant by virtue of the Maxwell equations dr{\/—gF^^) = 
dx{y/—gF^^) = 0. Again, the gauge held can only be a function of x and r. The boundary 
current, linearized about the gauge held and metric perturbations, is 


5F { 

a{g-x^F) r ~ ^ ^ 


V{dxH)5gtr + ( ^ 
ro 


VdrH 


59 


tx 


(5.31) 


As we are working at hnite temperature, there is another conserved quantity, Q, associated 
the the heat current. This conserved quantity is associated with a tensor ra 

G‘“' = V'T + “ 24,)FI “', (5.32) 

where is a Killing vector such that C^F = 0. Also, C^A = dijj and i^F = dcj), i^F being 
the interior product of ^ and F. With these dehnitions, it can be shown that 

= 3r , (5.33) 

provided that the Maxwell fl2.3p and traced Einstein fl2.2p equations are satished [T3] . 

An appropriate killing vector is = [1,0, 0,0], which satishes the re¬ 

quirement C^F = 0. There is a conserved quantity associated with fl5.32p as can be seen 
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by observing that is antisymmetric and that it is only a function x and r. Hence, 
dx{y/^G^^) = dr{y/^G^^) = 0, which follows from (15.331) . The conserved quantity is 
Q = yJ—gG™. As shown in [13], this quantity is the heat current in the boundary theory. 
Linearized about the perturbations, Q is 




5F 


2rL‘^ y a{ri — x^F) 


rL^a 




{d^Clx dxClr 


(5.34) 


+rL‘^a {drSgtx - d^Sgtr) + r^V[X q 1 


^^0 


drH + 2L‘^a - 


ro 


1 - 


^^0 


/ior^H 


H 


ro 


L^rdra + r^V HdrFl) 6gtx — ( L^rdxa + g,or^V ( 1 


^^0 


dxH + rV HdxH 6g, 


Hr 


In what follows, we will make use of the perturbative expansion in the disorder strength and, 

_2 

similarly to the initially uncharged case in section H] we will ignore terms of order V 5g and 
higher. Following m, the next step is to evaluate the constants fl5.3ip and Q fl5.34p near 
the event horizon at a hxed temperature T. Using the expansion for the metric coefficients 
fl5.9p - fl5.14p and the perturbation in fl5.17p - fl5.19p and fl5.22p - fl5.28p . we hnd 


r 


' S{o)F(^o) 

«(0)h(0) 


0(0) {E + dxar(o)) + ^ 



^9tx(o) 1 


(5.35) 


and 


Q = —2'kTLS gtx{o) = constant, 


(5.36) 


meaning that 5gtx{o) = constant. Next, we expand Q fl5.34p to next order in the near horizon 
expansion and use the linearized equations of motion to express it entirely in terms of near 
horizon data. This expression may then be used to find an expression for 5gtx{o) in terms of 
E. The next order expansion of Q is messy and contained in appendix]^ The result is that 


“(0) (fio - roVHt^i)) [E + dxar(o)} - ^TiTrodx (a(o)Sgtr(o)) 

, , o^0t7 zj ^0/^0 , 87r^T^roa(o)X(o) ^ . 

-I- I 87rTQ!(o) -l- 2—U/io77(i)-— I- -r2(o)(x) j i5gtx(o) — 0 , 


(5.37) 


where 11(0) (a;) is a function of metric coefficients near the horizon. It is found in appendix lAl 
To get an expression for in terms of E and 5gtx{o)^ we integrate fl5.36p over the common 
periodicity 2nN/kQ. Then, doing exactly the same thing with fl5.37p . we can relate E to 
^9tx{o)- Substituting this into the expression for gives us the desired relation between 
and E. The final result is 


r _ h 
'E~ h + h' 


(5.38) 


where 


ko pro/ig 

2tiN Jo [ 


o nr Stt^TVo 2 , n ^ a 

—- SnTa^o-) -——a(o)X(o) + ^(o)(a;) 


dx , 
(5.39) 
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h = 


X 


7 2 i'2-KN/ko 

^0 / 


4vr2iV2 j 
^27rA^/fco 


m 


«(0)<5(0)-^{0) 


(ix 


(5.40) 


2ro/ig 4rg/io 


L2 


L2 


VH{i) — StiTo 


( 0 ) 


Sti^T^vo 

m 


“(o)X(o) + ^(o)(a^) 


cix, 


U2 i‘2nN/ko _ 

^3 = 4^25^2JV2 Ko) “ ^0). 




( 5 . 41 ) 


All of the functions in the integrals fl5.39p - fl5.4ip may be expanded to second order in the 
disorder strength V, provided the solutions are known. The final step after the expansion 
is to take the disorder average fl2.8l) of fl5.38l) . The disorder averaged conductivity is then 
16 TrG]\f{a)£) = {J^/E)d. 


6 Summary and outlook 

In section |2] we construct a bottom-up holographic model with perturbatively charged dis¬ 
order. Starting from a bulk Einstein-Maxwell action, we include disorder in the dual theory 
by using a spectral representation. This technique is known to simulate a stochastic process 
[48] . This is achieved by including a randomly varying chemical potential (12.6p made up of a 
sum of N periodic functions along one of the boundary directions; an approach reminiscent 
of the disordered holographic superconductors studied in [17], [49] and [50] . The random 
chemical potential contains two parameters, a wavenumber fco which is held hxed in the large 
N limit and V which controls the strength of the disorder. The parameter V is taken to be 
small and is used as a perturbative handle to construct a bulk solution. A bulk gauge held 
is turned on which approaches the huctuating chemical potential fl2.6p near the spacetime 
boundary. By letting the gauge held backreact on the initially clean (i.e. zero disorder) AdS^ 
geometry, we construct an asymptotically AdS solution to the bulk equations of motion fl2.2p 
and fl2.3l) at second order in the disorder strength in section 12.11 

We evaluate the disorder average (12.8p of the second order metric coefficients fl2.18p . 
fl2.19p and fl2.20p . in section 12.21 and hnd that they are compactly expressed in terms of 
special functions. By carefully evaluating the large N limit, we hnd a divergence in the deep 
interior as r —)■ 00 in two of the metric coefficients fl2.2ip and fl2.22p . indicating that the 
solutions must be regulated. 

In section 131 we resum the disordered solution found in section [2] by adapting the standard 
Poincare-Lindstedt method for regulating perturbative solutions to diherential equations, 
similarly to the case of scalar sourced disorder in [52] • We hnd a regulated, second order, 
solution to the equations of motion with metric functions fl3.22p . fl3.23p and fl2.23p . where 
the previous noted divergences are removed. The averaged resummed solution is devoid of 
curvature singularities through second order in the disorder strength. 

With the resummed solution at hand, we study the resulting DC conductivity along the 
disordered direction in section [41 We directly access the DC conductivity of the model 
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by adapting a technique first proposed in PI- The idea is to turn on a source linear 
in time and take advantage of the existence of conserved quantities in the bulk to hnd a 
relationship between the boundary current and the magnitude of the applied electric held. 
The DC conductivity may then be extracted from Ohm’s law. By virtue of the nature of 
the disordered spacetime solution in this case, turning on a bulk gauge held perturbation 
along the disordered direction results in a complicated set of equations of motion. As is 
the case for single holographic lattices [B], the gauge held perturbation further sources a 
whole set of possible metric huctuations. By taking advantage of our perturbative handle, 
namely the disorder strength, the situation simplihes and the relevant metric huctuation can 
be solved via the linearized equations of motion up to second order in perturbations. The 
disorder averaged DC conductivity is computed in 04.111) which is found to be the usual AdS 
conductivity plus a correction at second order in the disorder strength, a result reminiscent 
of the single holographic lattices studied in [51]. 

Section [5] makes some observations about adding a disordered chemical potential to an 
initially clean system with a hnite charge density. In this case, the baseline geometry is 
the Reissner-Nordstrom-AdS (RN-AdS) solution. The presence of an initial charge density 
changes the character of the backreacted equations of motion, as the perturbative contribu¬ 
tions originating from the disorder source now mix with the baseline gauge held. The result 
is that geometry may receive corrections at all orders in the disorder source, unlike the ini¬ 
tially uncharged case studied in section [2l At second order, the traced Einstein equations are 
complicated and it is not clear that there is a compact analytic solution. It is still possible 
to extract some information about the form of the disorder averaged DC conductivity in this 
case by applying the techniques of [ 13 ] and HU The procedure is similar to the initially 
uncharged case in section 01 except that now the baseline solution (and hence the disordered 
solution) has an event horizon. By writing down a broad ansatz for a 3 -|- 1 dimensional dis¬ 
ordered geometry at hnite charge density and temperature (15.71) . we show that the disorder 
averaged DC conductivity may be expressed entirely in terms of near horizon data fl5.38|) . as 
is the case for the holographic lattices studied in HU- We leave the difficult task of hnding 
a disordered, hnite charge density, spacetime solution for future work. 

There many open questions with regard to explicit implementations of holographic disor¬ 
der, both in terms of studying the properties of the dual held theory as well as understanding 
the kinds of bulk geometries that arise in the process. 

A natural extension to our work here is to disordered holographic superconductors, such 
as those studied in m, Bsi and [50]. A spectral representation for the disordered chemical 
potential is also used in these studies and the properties of the superconducting transition are 
studied numerically. It would be interesting to understand how the inclusion of backreaction 
of the disorder source onto the spacetime geometry changes the picture here. For example, 
how is the appearance of the superconducting phase transition ahected? Is the critical tem¬ 
perature signihcantly changed? Is possible to get an analytical handle on a a critical amount 
of holographic disorder beyond which the conductivity becomes completely suppressed and 
the superconducting phase transition does not occur? This previous question is particularly 
pertinent with regards to many-body localization. If such a transition does occur, what kind 
of bulk geometry is required and how does it £t into conventional gravitational models? In 
particular, if a localization transition does occur in the dual theory, would the bulk probe 
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effectively become stuclfl at some radial position? Would this translate to a well defined mo¬ 
bility edge at the corresponding energy scale in the dual theory? It would also be interesting 
to understand the behaviour of time dependent probes in such a background. 

To fully study this problem, it may be necessary to move beyond the perturbative disorder 
studied in this paper and an analytical approach may be ill-suited and numerical solutions 
may be required. In such a case, it would be interesting to understand if the techniques used 
in [52] for scalar disorder and in for holographic lattices would be useful. 

In [52| , disorder is sourced by a scalar field in 2 -|- 1 dimensions and the disorder averaged 
metric is found to display an emergent Lifshitz scaling. It would interesting to classify 
the possible IR geometries that can be produced in this way, hopefully leading to a better 
understanding of disorder fixed points of condensed matter systems. For example, could an 
interior geometry with an emergent hyperscaling violation exponent be generated via a back- 
reacted disordered source? How about IR geometries that break rotations, i.e. in relation 
to Bianchi models [23], [21]? In order words, what kind of IR disordered hxed points can be 
constructed via holography? 

In [l6] , an ansatz for a disordered geometry is proposed and the dynamics of a scalar held 
in this background are studied. Using techniques from random matrix theory, a transition is 
observed which is reminiscent of a disorder driven metal-insulator transition. Understanding 
more than one explicit example of a geometry which displays this behaviour as well as the 
matter content required to support such solutions in a gravitational theory may shed light 
on the minimal ingredients necessary for accessing disordered phenomena via holography. 
It would also be worthwhile understanding how the proposed spacetime in [15] hts in with 
backreacted disordered geometries. 

Finally, our results for the initially uncharged case in sections [21l3]and|l]depend sensitively 
on the disorder distribution. In this paper we have focused on the effect of Gaussian random 
disorder. It would be interesting to extend our results to other distributions and understand 
how the resulting backreacted geometry and transport properties are modihed. 
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A Appendix 

As mentioned in section [5l one of the ingredients needed to compute the DC conductivity is 
the expansion of the heat current Q to order r — tq near the horizon tq. This results in the 
following condition 


a(o) (ho - roVHi^i)] [E + dxar(o)) - AirTrod^ (a(o)<5fi'tr(o)) 


(A.l) 


MTa{Q) + 


2rlnoV 

L2 




( 1 ) 


/igro , Svr^T^ro 2 
“(o)X(o) 


L2 


m 


^9tx(o) 


24ro7r^T^a(i) + 87r^TVoa(o)^ + 27rTroa{o)T^27rTroa(o) — 

^(0) ^(0) Vio). 


^9tx(o) 0. 


Using the linearized equations of motion, it is possible to solve for last term in brackets 
containing a(i), F(i), 5(i), p(i) and express it entirely in terms of horizon data with the result 
fl5.37p mentioned in the main text. To the relevant order in perturbations, we get 


D(o)(x) 


- 87iTL^a(^o)V{o) - 2/iol^?^oh(o)^^(i) + 327r2TVo?’o^a(o)^{o)h(o)^(o)^^(i) 

- 167r^T62/igroa(o)F(o)?7(o)5(o) + 1287r'^T^roFj)«?o)'^(o)X?o) “ 87r^T2LVo^(o)a(o)X(o)] 

+ 47rTroT^F(o)5(o) [l67r^T^F(o)Q!^o)'^^o)X(o)^xh(o) + 167r^T^F(o)Q!(o)?7(o)<^(o)X(o)^xa(o) 

- 327r^T2F(o)a^o)h(o)<^^o)^xX(o) + 167r^T2a^o)^(o)h(o)<5(o)'9x(5(o) + 2a(o)r7(o)^(o)4X(o) 

- «(0)^(0)X(0) <9x^(0) — ^V{0)^(0)X{0)dxOi(0) — «(0)h{0)X{0)^a:^(0)] 

+ levr^T^LVo [F(o)?7(o)(5fo)(^^«(o))^ + 9^(o)a?o)h(o)5(o)(45(o))^ + a^o)'^?o)(^xh(o))(45(o)) 

+ A'(0)«(0)<^fo)(^3;'^(0))(^xh(0))] — L‘^f'0 [h(0)^(0)(^x«(0))(^x^(0)) + «(0)^(0) (<9x^(0) )(^x^(0)) 

+'^(o)h(o)(^x^(o))^] “ 2L^roa(o)h(o)<^(o) [l67r^T^F(o)5^o)^x'^(o) + 167r^T^F(o)a(o)^(o)'9^<^(o) 
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